Molecule Formation in Optical Lattice Wells by Resonantly Modulated Magnetic 

Fields 
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We present a theoretical model for formation of molecules in an optical lattice well where a 
resonant coupling of atomic and molecular states is provided by small oscillations of a magnetic field 
in the vicinity of a Feshbach resonance. As opposed to an adiabatic sweep over the full resonance, 
this provides a coherent coupling with a frequency that can be tuned to meet resonance conditions 
in the system. The effective Rabi frequencies for this coupling are calculated and simulations show 
perfect Rabi oscillations. Robust production of molecules with an adiabatic sweep of the modulation 
frequency is demonstrated. For very large oscillation amplitudes, the Rabi oscillations are distorted 
but still effective and fast association is possible. 

PACS numbers: 03.75.Nt, 03.75.Ss, 36.90.+f 



Creation of molecules from ultracold atoms has been 
successfully achieved by making an adiabatic sweep of the 
magnetic field across a Feshbach resonance UllZklQLbl 
6|. Feshbach resonances can be used to create molecules, 
because of the transformation of a scattering state into 
a new bound state, when the scattering length changes 
sign. This occurs when the magnetic field varies near a 
Feshbach resonance 71: 



a sc (B) = dbg ( 1 



B — Bo 



(1) 



Here ab g is the background value of the scattering length, 
and Bo and A are the location and width of the reso- 
nance. In a recent publication y| it was demonstrated 
that near Feshbach resonances, molecules can be created 
by harmonic modulation of the magnetic field: 



B(t) =B' + bsm(u) B t), 



(2) 



where is resonant with the energy difference between 
the atomic and the molecular state at the magnetic field 
B'. It is the purpose of this article to analyze this pro- 
cess for two atoms in an optical lattice well, which is 
a very clean system for which many-body effects can 
be neglected. One argument in favour of producing 
the molecules in the lattice resonantly rather than by 
sweeping across the Feshbach resonance is that it enables 
the preparation of well-controlled coherent superposition 
states of atomic and molecular components. 

Ultracold molecules in optical lattices have been cre- 
ated by photoassociation 0, , and proposals for cre- 
ating heteronuclear molecules have also been presented 
|H|. In [13 , Rabi oscillations have been seen with this 
method. In |12.ll3| the production of molecules via Fes- 
hbach resonances in a ID confining lattice potential has 
been studied and very recently the first successful cre- 
ation of molecules in separate wells of a 3-dimcnsional 
optical lattice by sweeping across a Feshbach resonance 
has been reported 01- As ^ turns out these molecules 



can be very long-lived (up to 700 ms) since they do not 
suffer from decay due to inelastic collisions. 
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Figure 1: The s-wave spectrum of the potential around the 
Bo = 155 G Feshbach resonance in 85 Rb in a harmonic oscilla- 
tor trap with u — 2n- 30 kHz. The variation of the scattering 
length a sc as a function of B is shown as a dot-dashed line. 
Note that the bound (molecular) state with negative energy 
occurs in the a > (repulsive effective interaction) regime. 
The typical adiabatic production of molecules by a magnetic 
field scan across the resonance is indicated by the arrow (1). 
Arrow (2) illustrates the production of molecules by resonant 
association as analyzed in the present article, while the dotted 
arrows (3) show an accidental three-level resonance (discussed 
in the text). 



In [3] a commonly used resonance in 85 Rb is applied. It 
has the parameters Bo=155.041 G, A =10.71 G, and 
abg = —443 ao, where ao is the Bohr radius. This reso- 
nance has previously been used to produce molecules 
and to tune the scattering length to make stable atomic 
BECs of 85 Rb atoms 

We now look at a simple two-atom model for molecule 
formation in a three-dimensional optical lattice which is 
sufficiently deep that the atoms do not tunnel between 
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lattice sites. The system may either be in the Mott insu- 
lating state [l7| , or we may have a mean field state with 
number fluctuations in the site occupancy, in which case 
the collisional interactions apply to the two-atom compo- 
nent of the state on each site [l8j . The assumption is that 
we have two identical atoms in a lattice well and neglect 
couplings to atoms in other wells. The lattice potential 

is V(x,y,z) = VoJ2i=x> y ,z sin2 ( kx i)> k = X> where 
A is the wavelength of the standing wave laser field of 
the lattice, and we make the harmonic and isotropic ap- 



proximation in each well: V(x, y, z) 



The potential for two atoms in this harmonic potential 
decouples into parts describing the center-of-mass and 
the relative motion. The interaction between atoms is 
modelled by the Huang s-wave pseudopotential and 
the total potential for the relative motion becomes 



V(r) = i(m/2)wV 



4irh 2 a sc (B) 



(3) 



This is an adequate approximation to the more compre- 
hensive analysis [2(j which takes into account the spatial 
confinement and the energy dependence of the scattering 
length. We note that a more precise treatment amounts 
in principle to the application of the exact energy eigen- 
values and matrix elements instead of the analytical ex- 
pressions offered by the pseudo-potential in J3J). Such an 
analysis will lead to the same physical behaviour of the 
system. For our examples with Rb atoms and typical op- 
tical lattice parameters, the characteristic values of the 
scattering length is about one order of magnitude smaller 
than the harmonic oscillator length, and we expect that 
the analytical model provides a reasonable approxima- 
tion to the exact dynamics |2*l) . 

We adopt harmonic oscillator units y/H/ (m/2)u> for 
length and fko for energy, and we express the energy of 
the s-wave states by a generalized harmonic oscillator 
quantum number: E t (B) = 2v l (B) + 3/2. Following [H 
|23| . Vi(B) are obtained as the solutions of the equation 



2r(-i/<) 



Y{-Vi - 1/2) a sc (B) 



(4) 



where T is the gamma function. The energy spectrum is 
plotted in Fig. the 85 Rb resonance. 

Assuming a time dependent magnetic field, the eigen- 
states depend on a sc and thus on t. We shall label the 
states \ip n {t) >, n — 0,1,2,... with increasing energy. 
They will at every instant form a complete set. Expand- 
ing the wave function \ip(t) > on these states we obtain 
the following equations for the expansion coefficients c m : 



ih dCrJA = Em{t)Cm{t) _ th f- dO mn (t,t>) 



df 



dt' 



Cn(t) 

(5) 



where O mn (t, t') =< il} m {t)\ip n (t!) >• Introducing the 
v- values v m and v n corresponding to ij} m {t) and tp n (t') 



respectively and by following a trick in |22j , or by explic- 
itly evaluating the integral of the wave functions [24[ , we 
obtain 



T(-v m )Y(-v n - i) - r(-f n )r(-z/ m - 1) 



\Jv{-v m - \)T{-v m ) ■ yjv{-v n - \)T{-v n ) 
f{vm)f{y n ) ( a sc {t') a sc (t)\ 



where 



a sc {t) V a sc (t') 



*(-!/)-*(-!/- |) 



(6) 



(7) 



and ^{v) = ■4-\n{r{v)) is the digamma function. 

Inserting now the oscillating magnetic field J3J in J5J 
we get the following coupling between the states: 

dc (t) °° 
ih — ^ — = E m (t)c m (t) + ihcos(u} B t) V" U, mn {t)c n {t) 



n=0 



(8) 



where the Rabi frequencies O mn are given by 



-(6/A)c B A 2 2 dO„ n 

TJV'n) n 



(B-B a )(B-B - A)' 



(9) 



The Rabi frequencies are proportional to 6/ A, and the 
rest of the expression is plotted in Fig. [21 We remark 
that Q nn = 0. Equivalent to the radiative coupling of 
atomic energy levels, if Hlob matches the energy difference 
between two specific levels m and n, Rabi oscillations 
between these two levels with a constant frequency VL mn 
are expected as long as & C A. Notice that due to the 
factor ws in Eq. El there is a much stronger resonant 
coupling to the lowest (molecular) state than among the 
atomic states on the molecular side of the resonance. 

Using Eq. JSJ we have calculated the dynamics of the 
resonant association of two 85 Rb atoms in an uj — 2tt ■ 30 
kHz lattice well. Referring to the above resonance we use 
B' = B Q + 0.4A = 159.33 G, and b = 0.01 and we start 
out in the lowest atomic state (on the 2 nd lowest curve 
in Fig. The expected resonant frequency is found 

from Q: u™ s = (E-l(B') - E (B'))/h = 5.079 u = 2n ■ 
152.4 kHz. Although the calculation is done with 10 basis 
states the dynamics in Fig. |3 shows full two-level Rabi- 
oscillations as one would expect. A fit of the molecular 
state population to a sin 2 (f2/2 • t) function shows that 
the Rabi frequency is = 2tt ■ 1.84 kHz which exactly 
matches the value of the Rabi frequency displayed in Fig. 
121 Matching the oscillation time such that £lt — it one 
can transform the atomic pair into a molecule with unit 
probability, and using the interaction time as a control 
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Figure 2: Amplitude normalized Rabi frequencies on reso- 
nance for 85 Rb. The parameters are as in Fig. Q The states 
are labelled 0, 1, 2, . . . with increasing energy. E.g., fiio is the 
coupling between the two lowest states. 




Figure 3: Dynamics of the resonant process at B' = Bo + 
0.4A. The field is modulated at an angular frequency ujb — 
2iv ■ 152.4 kHz, with the amplitude b = 0.01A. The other 
parameters are as in Fig. 



parameter, a wide range of well-controlled superposition 
states of the molecular and atomic component can also 
be prepared. 

At certain magnetic fields three levels are coupled in 
a ladder configuration, because hojg s matches not only 
the energy distance to the molecular state, but also the 
distance to a higher lying atomic state (see Fig. This 
leads to a more complicated three-level dynamics, but 
since the coupling to the higher lying atomic state is 
much smaller than the coupling to the molecular state 
(see Fig. the higher lying atomic level becomes only 
sparsely populated. This is fortunate from the point of 
view of molecular formation. The lowest B-field on the 
molecular side of the resonance where such a three-level 
resonance occurs is B' = Bo + 0.3582A = 158.88 G shown 
with the dashed arrows in Fig. n(c^ s = 2tt ■ 125.1 kHz). 
Simulation of the dynamics in this case shows oscillations 
with a peak population of 6 % in the excited atomic 



state. All remaining population is distributed between 
the molecular and the lowest atomic state. At the next 
three-level resonance at B' = B Q + 0.4385A = 159.74 G 
the peak population of the higher lying atomic state be- 
comes only 4 %. 

The scheme should work just as well for Feshbach res- 
onances in other species like 87 Rb although the Feshbach 
resonances are more narrow and situated at higher mag- 
netic fields. According to [25| the broadest resonance in 
87 Rb occurs between atoms in the F=l mp=l hyperfine 
state and has the following parameters: Bq =1007.60(3) 
G, A =0.20(3) G, cibg = 100. 5a . This resonance has 
been used to produce molecules in an optical dipole trap 
@ and it was also used in the recent optical lattice ex- 
periment |l4j . Some Rabi frequencies for this system are 
shown in Fig. The couplings are of the same order of 
magnitude, but the plot is qualitatively different from the 
plot in Fig. [21 The difference can be partly attributed 
to the different behaviour of the scattering length as B 
approaches Bq — A (positive scattering length side in 
87 Rb) compared to when B approaches Bq + A (positive 
scattering length side in 85 Rb) - see QJ. In the latter 
case the scattering length approaches and the resonant 
frequency ujg s entering diverges for coupling to the 
molecular state, because the binding energy of the molec- 
ular state in the pseudopotential model varies as 
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2(m/2)al ma* g (l - A/(B ~ Bo))* 
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Figure 4: Amplitude normalized Rabi frequencies in 87 Rb (cf. 
Fig. Here the positive scattering length (molecular) side 
is to the left of the resonance, because ab g = 100.5 ao > 0. 
We assume the same trapping frequency u — 2tt- 30 kHz as 
above. 



Application of a 7r-pulse in the above scheme relies on 
the resonance and Rabi frequencies being exactly known 
and identical for all lattice sites. A more robust scheme 
is to sweep the frequency slowly across the resonance 
frequency. Using a constant amplitude linear frequency 
sweep from 10% below to 10% above cj^ s , we get the 
dynamics shown in Fig. [S] After just a few Rabi oscilla- 
tion periods the molecular state population settles near 
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100 %. In addition to correcting for inhomogeneity ef- 
fects, a frequency sweep will also introduce some robust- 
ness against technical noise in the magnetic field value, 
as long as this noise is present only on longer timescales 
than that given by the sweeping time. This is especially 
interesting regarding applications to more narrow Fesh- 
bach resonances like in 87 Rb. 
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Figure 5: 5 ms linear frequency sweep from 10% below to 
10% above the resonance frequency ufg s — 2n ■ 152.4 kHz at 
an applied field B' — Bo + 0.4A. The modulation has the 
constant amplitude: b = 0.01A. 



One might also consider the possibility of using larger 
oscillation amplitudes. This makes the association 
quicker, can be used to keep the magnitude of the har- 
monic modulation well above the technical noise level 
and gives also some robustness against a noisy value of 
B' . Our simulations show that if the coupling strength is 
not much smaller than A, the resonance frequency is sig- 
nificantly shifted and the oscillations are distorted, but 
they show a surprisingly high molecular peak population. 
For instance with B' = Bq + 0.4A and b = 0.2A we must 
use a value of lob which is 45 % greater than that in Fig. 
|3]to get maximum population oscillations, and we then 
get a molecular peak population of around 80 % (Fig. 
EJ). We attribute the frequency shift to the strongly in- 
creasing slope of the molecular level curve (Fig. as the 
B-field is modulated to larger values, farther away from 
the resonance. 

To summarize, we have analyzed the formation of 
molecules from atoms in an optical lattice by a resonant 
modulation of a B-field which causes Rabi-oscillations 
between the atomic and the molecular state in the vicin- 
ity of a Feshbach resonance. For identical atoms, the 
problem separates in a center-of-mass and a relative mo- 
tion part, and the relevant s-wave scattering part was 
solved exactly in terms of known analytical solutions of 
the time-independent Schrodinger equation. Complete 
conversion to molecules was seen both with a 7r-pulse of 
the interaction and with a frequency swept interaction. 
In the case of different atoms, relevant for the formation 
of heteronuclear molecules, we expect the coupling be- 
tween center-of-mass and relative motion to cause highly 
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Figure 6: Dynamics of the process with a very high oscilla- 
tion amplitude: b = 0.2A. B' = B + 0.4A. The resonance 
frequency has shifted 45 % compared to the value in Fig. |3]so 
ujb = 1.45 ■ 2ir ■ 152.4 kHz. The frequency of the population 
oscillation is 33 kHz whereas the Rabi frequency predicted 
from © is 37 kHz. There are 20 basis states included in the 
calculation. 



non-trivial effects. 

In the experiment the large number of 85 Rb atoms 
is not pairwise confined in a lattice potential, and 
rather than perfect Rabi oscillations between atoms and 
molecules, the experiments show small rapid oscillations, 
dying out towards a steady state with about 30 % of 
the atoms converted into molecules. Many-body effects, 
including collisions in the gas, probably cause the dis- 
crepancy with the simple two-atom model. The simple 
picture remains, however, useful for the identification of 
the resonant process in the system, and we imagine that 
the swept modulation frequency, proposed in this article, 
may lead to enhanced molecular fractions in the many- 
body system. 

Since efficient conversion into molecules inside an op- 
tical lattice was recently demonstrated by an adiabatic 
sweep of the B-field across the full 87 Rb Feshbach reso- 
nance |14| , let us conclude by pointing out some relevant 
features of the modulated field approach. Our approach 
is equivalent to driven two-level transitions in laser and 
NMR spectroscopy, and by varying the phase, amplitude 
and duration of the applied B-field modulation, we have 
more degrees of freedom to control the transition process. 
For a number of quantum optics studies, e.g., for quan- 
tum non-demolition photon counting and non-classical 
state preparation and for quantum state tomography, the 
coherent Rabi oscillation phenomenon is extremely use- 
ful, and analogous proposals for trapped atoms [2|| can 
be implemented with the modulated B-fields. In con- 
nection with applications to quantum information, uni- 
tary operations are requested that can be applied to un- 
known initial quantum states, and for this purpose single 
pulses, or possibly composite pulse sequences with built- 
in robustness against frequency and amplitude errors, are 
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more versatile than the adiabatic passage across the full lai Nygaard, and Uffe Vestergaard Poulsen for comments 
resonance. on the manuscript. 
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